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Parameters  representing  physical  measurements  Involved  In  the  dissolved 
oxygen  deficit  Equation  (1) : 

V - Plow  rate  of  the  river  before  entering  Reach  1,  In  10 
*0  gallons/day 

L.  - Initial  BOD  level  of  the  stream  before  entering  Reach  1, 

“O  In  mg/1 

D.  - Initial  oxygen  deficit  of  the  stream  before  entering  Reach  1, 
**0  in  mg/A 

K - Deoxygenation  constant.  In  day  ^ 
r - Reaeration  constant,  in  day  ^ 
t - Flow  time  along  the  reach.  In  days 
F - Fraction  of  river  bottom  covered  with  sludge 
Q - A coefficient  in  the  sludge  term,  determined  empirically 
In  day  ^ 

K - Oxygen  uptake  rate  per  unit  area  of  stream  bottom  surface, 

® 2 

in  gm/m  /day 

R^  - Hydraulic  radius  of  the  river  cross  section,  in  meters 

V - Volume  of  the  effluent  released  Into  the  river.  In  10^ 

gallons/day 

E - BOD  concentration  of  effluent,  before  treatment.  In  mg/1 

Aside  from  the  above  parameters  which  are  Involved  in  the  oxygen 
sag  equation  of  each  reach,  the  following  parameters  will  also  enter  Into 
the  formulation  of  the  model  at  various  stages. 

Parameters  representing  quantities  determined  by  management  decisions  or 
feasibility  considerations: 

S - Maximum  allowable  oxygen  deficit.  In  mg/1 

jp  - Minimum  required  or  feasible  BOD  removal  by  a specified 

sequence  of  treatment  components  In  a given  reach,  a fraction 

T - Maximum  possible  or  feasible  BOD  removal  by  a specified 

sequence  of  treatment  components  in  a given  reach,  a fraction 

pp  - Minimum  required  or  feasible  BOD  removal  by  a specified 
single  treatment  component  In  a given  reach,  a fraction 
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The  existing  approaches  to  the  mathematical  modeling  and  optimization 
of  a water  pollution  problem  are  briefly  surveyed.  A proposed  geometric 
programming  model  of  a water  pollution  and  treatment  system,  which  easily 
lends  Itself  to  the  theory  and  application  of  nonlinear  programming  sensi- 
tivity analysis  techniques,  is  studied  in  detail.  As  in  previous  work  by 
the  author  of  the  model,  the  optimal  waste  treatment  facilities  along  the 
Upper  Hudson  River  are  presented  for  three  different  environmental  policies, 
leaving  the  sensitivity  analysis  study  for  a following  report.  Previous 
results  were  obtained  using  a geometric  programming  code,  while  the  present 
study  makes  use  of  the  SUMT  code.  It  is  shown  that  the  policy  of  fixed 
dissolved  oxygen  requirement  yields  the  minimum  annual  waste  treatment  cost, 
relative  to  two  other  policies.  Moreover,  it  is  shown  that  a variable 
dissolved  oxygen  policy  yields  relatively  uniform  treatment  levels  in  the 
treatment  plants  at  appreciably  reduced  cost  compared  to  the  costs  Involved 
in  strict  uniform  treatment  policy.  The  results  of  the  present  study  are 
consistent  with,  but  not  identical  to,  the  findings  reported  in  previous 
work.  The  discrepancy  in  the  coefficients  involved  in  the  dissolved  oxygen 
constraints  used  by  these  two  studies  is  the  main  reason  for  the  differ- 
ences observed.  A listing  of  a computer  program  developed  to  calculate 
these  coefficients  is  included. 
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20.  Abstract  - continued 

Reaving  the  sensitivity  analysis  study  for  a following  report. 

Previous  results  were  obtained  using  a geometric  progranmlng  code,  while 
the  present  study  makes  use  of  the  SUMT  code.  It  Is  shown  that  the 
policy  of  fixed  dissolved  oxygen  requirement  yields  the  minimum  annual 
waste  treatment  cost,  relative  to  two  other  policies.  Moreover,  It  is 
shown  that  a variable  dissolved  oxygen  policy  yields  relatively  uniform 
treatment  levels  In  the  treatment  plants  at  appreciably  reduced  cost  com- 
pared to  the  costs  Involved  In  strict  uniform  treatment  policy.  The 
results  of  the  present  study  are  consistent  with,  but  not  identical  to, 
the  findings  reported  in  previous  work.  The  discrepancy  in  the  coeffi- 
cients involved  in  the  dissolved  oxygen  constraints  used  by  these  two 
studies  is  the  main  reason  for  the  differences  observed.  A listing  of 
a computer  program  developed  to  calculate  these  coefficients  is 
included. 
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1.  Introduction 

The  present  paper  reports  on  the  solution  and  analysis  of  a geometric 
programming  water  pollution  control  model  [1]  formulated  by  J.  G.  Ecker  and 
applied  by  him  to  a problem  involving  data  for  the  Upper  Hudson  River.  A 
basic  premise  motivating  the  environmental  pollution  problem  is  that  the 
pollutant  level  should  not  "significantly"  diminish  the  "quality  of  life." 
Stated  more  simply  and  affirmatively  in  the  context  of  an  important  criterion 
of  stream  water  pollution,  the  dissolved  oxygen  of  the  stream  should  be  kept 
above  a certain  level.  The  standard  for  an  acceptable  level  of  dissolved 
oxygen  in  a stream  is  usually  set  by  the  various  environmental  protection 
agencies.  It  is  important  to  note  that,  other  than  dissolved  oxygen  (DO), 
there  are  many  other  possible  measures  of  water  quality,  e.g.,  temperature, 
caste,  color,  odor,  dissolved  chloride,  or  turbidity.  However,  the  dissolved 
oxygen  appears  to  be  the  most  common  water  quality  measure.  The  reason  for 
this,  as  Sobel  [2]  points  out,  may  be  paraphrased  as  follows: 

Dissolved  oxygen  fluctuates  in  response  to  many  of  the 

phenomena  reflected  in  the  other  measures;  and  it  may  serve 
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as  a surrogate  for  them.  DO  Improvement  programs  often  have 
side  benefits  associated  with  the  Improvements  In  all  the 
measures  mentioned.  The  DO  content  of  the  water  plays  an 
Important  role  In  supporting  aquatic  life.  Its  control  may 
be  accomplished  In  a variety  of  ways  that  have  motivated 
rather  sophisticated  problem  formulations. 

A major  cause  of  the  violation  of  dissolved  oxygen  standards  In  streams 
Is  the  release  of  Industrial  wastes.  Those  wastes  usually  contain  organic 
materials  which.  In  the  course  of  decomposition,  consume  and  diminish  the  dis- 
solved oxygen  content  of  the  water  and  endanger  aquatic  life  and  the  quality 
of  the  surrounding  environment. 

A standard  measure  for  organic  waste  quality  Is  the  level  of  Biochemical 
Oxygen  Demand  (BOD)  concentration,  defined  as  the  quantity  of  dissolved  oxygen 
In  the  stream  required  to  stabilize  the  organic  material  In  the  discharge  over 

O 

a five-day  period  at  20  c. 

In  order  to  meet  a specified  dissolved  oxygen  level,  the  discharged 
waste  volume  and  Its  BOD  concentration  must  obviously  be  limited.  Fortunately, 
the  organic  wastes  released  Into  the  streams  do  not  consume  oxygen  indefinitely. 
This  Is  due  to  the  fact  that  these  materials  eventually  decompose  and  stabilize. 
Moreover,  due  to  reaeratlon  along  the  stream,  the  consumed  oxygen  is  constantly 
replenished  by  absorption  of  oxygen  from  the  air. 

Using  the  dissolved  oxygen  content  as  a pollution  criterion,  numerous 
attempts  have  been  made  to  model  and  control  the  quality  of  a stream,  using 
mathematical  programming.  Such  attempts  have  generally  strived  for  one  or 
both  of  the  following  goals: 

1)  Provide  guidelines  for  allocation  of  treatment  facilities 
at  each  discharge  point. 

11)  Determine  the  optimal  waste  treatment  level  at  each  discharge 
point . 

These  gosls  have  to  be  met  such  that  a certain  dissolved  oxygen  standard  Is 
achieved  with  a minimal  expenditure  of  resources. 

It  Is  clear  that  the  DO  profile  equation  along  a stream  and  the  waste 
treatment  cost  function  play  important  roles  In  any  water  quality  management 
model.  The  nonlinearity  of  the  dissolved  oxygen  profile  (oxygen  sag) 
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Equation  [3]  and  treatment  cost  function  [4],  [5j,  [6]  render  such  models 
Inherently  nonlinear.  However,  In  the  early  attempts  by  Delnlnger  [7]  and 
Delnlnger  and  Sobel  [8],  the  problem  was  cast  as  a linear  programming  model. 
Delnlnger  structured  linear  programming  models  which  utilize  various  approxi- 
mations of  the  differential  equation  used  to  describe  the  dissolved  oxygen 
profile  of  the  streams.  Following  such  pioneering  efforts,  the  problem  was 
further  Investigated  In  the  format  of  a linear  programming  model  by  Sobel  [2], 
Kerri  [9],  Revelle,  Loucks,  and  Lynn  [10],  [11],  Graves  and  Hatfield  [12], 
Arbabl  and  Elzlnga  [13]  and  others.  As  McNamara  [14]  pointed  out,  however, 
linear  programming  formulations  unduly  restrict  the  form  of  the  functional 
relationships  Involved  In  the  model. 

A large  number  of  nonlinear  models  were  subsequently  developed  to 
optimize  the  abundant  variations  of  the  problems  that  arise  In  water  quality 
management.  Llebman  and  Lynn  [15],  Shlh  and  Krlshnan  [16],  Keegan  and  Leeds 
[17],  Meier  and  Beightler  [18],  [19]  and  Dysart  [20]  mainly  applied  dynamic 
programming  to  model  the  problem. 

Ecker  and  McNamara  [21],  Ecker  [1]  and  McNamara  [14]  were  able  to 
formulate  an  Important  class  of  pollution  control  and  facility  design  problems 
as  geometric  programming  models.  In  [21]  they  use  the  data  provided  in  the 
dynamic  programming  approach  of  Shlh  and  Krlshnan  [16]  and  obtain  essentially 
the  same  optimal  design.  Among  many  appealing  features  of  the  Ecker  models 
are  that  they  readily  lend  themselves  to  the  sensitivity  analysis  study 
approach  motivated  by  the  work  of  Flacco  and  McCormick  [22]  and  developed  by 
Flacco  [23].  Ecker  [1]  points  out  that  "the  model  Is  particularly  useful 
for  sensitivity  analysis  Involving  changes  In  the  stream  standards." 

The  present  paper  reports  on  the  application  of  the  above  model  for 
the  optimal  facility  design  and  treatment  level  allocations  along  the  Upper 
Hudson  River,  for  three  alternative  pollution  control  policies.  As  noted, 
these  problems  have  been  solved  by  the  author  of  the  models  using  a geometric 
progranmlng  code,  while  the  present  study  makes  use  of  a general  nonlinear 
programming  formulation  and  code,  SUMT  [24].  A computer  routine  was  developed 
to  calculate  the  coefficients  Involved  In  the  oxygen  deficit  constraints 
and  the  resulting  data  manipulations  have  yielded  somewhat  different  (and 
apparently  more  accurate)  coefficients  than  those  used  In  the  previous  study. 
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The  results  of  the  present  study  show  that  for  the  Upper  Hudson  River 
the  optimal  annual  pollution  control  cost  Is  minimum  under  the  fixed  dis- 
solved oxygen  requirement  policy  and  maximum  under  the  uniform  treatment 
requirement  policy.  The  third  policy,  which  allows  a variable  dissolved 
oxygen  requirement,  yields  the  same  stream  dissolved  oxygen  concentration 
as  the  uniform  treatment  policy,  but  under  considerably  reduced  annual  treat- 
ment costs.  Moreover,  this  policy  yields  relatively  uniform  treatment  levels 
that  are  easier  to  administer  than  the  other  policies.  These  results  are 
consistent  with,  but  not  Identical  to,  the  findings  reported  In  [1]. 

2.  Model  Description 

2.1  Background . As  mentioned  earlier.  In  order  to  meet  the  dissolved 
oxygen  standard  set  by  environmental  rules,  the  discharged  waste  volume  and 
Its  BOD  concentration  Into  the  streams  must  be  limited.  To  maintain  a tolerable 
release  of  BOD  concentration,  an  Industry  must  first  process  wastes  In  treat- 
ment plants.  Figure  1 depicts  a typical  arrangement  of  the  facilities  In  a 
treatment  plant.  As  noted,  the  segment  of  the  stream  from  one  treatment 
plant  to  the  next  Is  called  a "reach."  (More  generally,  a reach  is  taken  to 
be  that  portion  of  a stream  along  which  there  Is  not  a sudden  and  significant 
change  In  the  water  quality  and  stream  parameters.)  This  allows  one  to  divide 
the  stream  under  study  Into  various  segments  (reaches)  and  use  fixed  parameters 
to  define  the  oxygen  profile  along  each  segment.  As  Indicated  In  Figure  1, 
the  Initial  treatment  process  Is  called  a "primary  treatment,"  Intermediate 
processes  are  called  "secondary  treatments"  and  the  last  process  Is  called 
a "final"  or  "advanced  treatment"  process.  Primary  and  secondary  treatments 
taken  together  are,  albeit  somewhat  ambiguously,  referred  to  as  "total 
secondary  treatments." 

Depending  on  the  BOD  discharge,  length,  and  other  characteristics  of 
the  reach,  the  minimum  oxygen  deficit  may  either  occur  at  the  end  of  each 
reach  or  within  the  reach,  as  shown  In  Figures  2a  and  2b,  respectively.  In 
the  formulation  of  the  model  under  study  It  Is  assumed  that  the  DO  profile 
takes  the  form  given  In  Figure  2a. 

To  meet  environmental  standards,  the  dissolved  oxygen  along  the 
stream  must  at  all  times  be  above  the  minimum  allowable  dissolved  oxygen 
level  dictated.  In  terms  of  oxygen  deficit,  the  oxygen  deficit  at  any 
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Figure  2a.  Dissolved  Oxygen  Profile  Along  Reach  i. 


Natural  Dissolved  Oxygen  Level 
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time  and  at  any  point  along  the  stream  should  not  exceed  the  maximum 
allowable  oxygen  deficit  dictated  by  the  environmental  standards. 

The  equation  for  the  dissolved  oxygen  deficit  at  time  t along 

a particular  reach  Is 

°t  “ *^1^  *^2°  *^3 

where 

L >■  BOD  concentration  of  the  stream  at  the  top  of  the  reach 
(after  waste  release) 

D - Dissolved  oxygen  deficit  at  the  top  of  the  reach,  and 
c^  , c^  and  c^  are  coefficients  which  depend  upon  time 

and  reach  parameters  to  be  discussed  in  the  next  section. 


The  BOD  concentration  L of  the  stream  at  the  top  of  each  reach  can 
be  calculated  from  the  following  simple  mass  balance  equation 


L 


V +V 
E R 


e + 


v„+v„ 

E R 


(2) 


where 


■ Volume  of  effluent  per  day  into  the  stream 

e - BOD  concentration  of  effluent  into  the  stream 

- Volume  of  stream  per  day  before  entering  the  reach 

* BOD  concentration  of  stream  at  the  bottom  of  the  reach 


The  BOD  concentration  at  the  bottom  of  Reach  1 can  be 

calculated  If  the  BOD  concentration  at  the  top  of  Reach  1 Is  known, 
as  follows 


h,.  - % "l 


(3) 


where  c^  Is  a coefficient  which  depends  upon  time  and  reach  parameter. 


Hence,  for  a sequence  of  reaches,  if  we  know  the  BOD  concentration 
and  oxygen  deficit  at  the  top  of  the  first  reach,  we  can  readily  obtain  the 
oxygen  deficit  profile  along  the  entire  stream  (for  all  reaches),  provided 
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that  reach  parameters,  stream  volume,  effluent  volume,  and  waste  discharged 
by  each  treatment  plant  are  known. 


2.2  The  Model.  The  model  under  study  is  a geometric  prograomlng 
model  developed  for  allocating  treatment  requirements  along  the  stream  so  as 
to  meet  dissolved  oxygen  standards  for  each  reach,  while  minimizing  the  total 
annual  cost  of  all  treatment  activities. 


The  problem  variable  t^^  Is  defined  as  the  fraction  of  BOD  remaining 


In  the  effluent  after  It  passes  through  the  jth  treatment  process  of  Reach  1. 
Figure  3 Is  a schematic  depiction  of  the  treatment  facilities  for  Reach  1 of 
the  Upper  Hudson  River. 


hi 


BOD  Concentration 


®1  hi 


h ^11  h2 


h ^11  *^12  *^13 
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is  the  BOD  concentration  of  the  waste  of  volume  v entering 

the  waste  treatment  system  1,  so  concentration  of  the 

treated  waste  remaining  after  process  1 and  entering  process  2, 

leaves  process  2 and  enters  process  3 and  ^]^i^X2’'13^1  process  3 

and  enters  the  stream.  As  mentioned  before  in  the  formulation  of  the 
model,  it  is  assumed  that  the  maximum  oxygen  deficit  within  each  reach  will 
occur  at  the  end  of  the  reach. 

It  is  stipulated  that  the  yearly  operating  cost  c 


ij 


of  process  J 


in  Beach  1 depends  upon  t^^  as  follows 

'ij  - 

where  c^^  and  a^j  sre  nonnegative  constants  that  can  be  determined  once 

the  kind  of  treatment  facility  is  specified.  If  there  are  n reaches  and 
m^  processes 
minimized  is 


m^  processes  in  reach  i then  the  total  yearly  operating  cost  F to  be 


n 1 


I I t i 

1-1  j-1 


-a 


ij 


(4) 


which  constitutes  the  objective  function  of  the  model. 

Theru  are  four  types  of  constraints  involved  in  the  formulation  of  the 


model: 


a.  Dissolved  oxygen  deficit  constraint. 

Dissolved  oxygen  deficit  along  teach  1 must  be  less 
than  maximum  allowable  oxygen  deficit  s^  along  Reach  1. 
This  constraint,  after  appropriate  manipulation,  is 

■l  *2  "l 

“11  'ij  “12  'll  * •••  * “11  'ij  i ' 

where  the  positive  constants  depending  on  s^ 

and  the  specific  parameters  of  Reach  1 and  those 
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upstream  of  Reach  1.  For  illustrative  purposes,  the 
dissolved  oxygen  deficit  constraints  of  Problem  for 


Reaches  1 and  2 and  the  calculation  of  the  u 
these  reaches  are  given  In  Appendix  1. 


IJ 


for 


b.  Combined  treatment  requirement. 

The  fraction  of  BOD  removed  by  a sequence  of  treatment 
components  In  Reach  1 may  be  bounded  from  above  or  below 
by  specified  fractions  P^  or  P^  respectively,  yielding 

constraints  of  the  form 


or 


1-  " t IP 
JcK* 


^ ^ ^li  — -1 

jeK* 


(6) 


I i 


where  K*  is  the  set  of  indices  of  the  subject  treatment 
components  In  Reach  1. 

c.  Operating  range  constraint  for  components  of 
treatment  facilities  In  Reach  1. 

The  fraction  of  BOD  removed  by  component  j in  Reach  1 
may  be  bounded  from  above  or  below  by  specified  fractions 


PP 

ij 

form 


or 


or  PPj^j  , respectively,  yielding  constraints  of  the 


1 - t < PP 
ij  - ij 


d.  Natural  constraints. 

The  variable  t^j  , the  fraction  of  BOD  remaining  in  the 

process  in  leach  1 after  process  j , is  a nonnegative  fraction, 
thus 
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0 ^ ^ 1 for  all  i,j  . 

After  some  simple  algebraic  manipulations,  the  given  model  for  a 
stream  Involving  n reaches,  with  Reach  1 containing  processing 

components,  can  be  formulated  as 

m. 

n i -a. 

minimize  F » ^ I c.,  t / 

1-1  j-1 


subject  to 

a.  Dissolved  oxygen  deficit  constraints  (5) 


“l 

n 

t,  ^ 

< 

1 

j-1 

ij 

“2 

"2 

— 

n 

t, , 

+ 

u„„ 

n t„,  < 

1 

J-1 

ij 

22 

• 

3-1  ' 

m. 

• 

• 

m 

1 

2 

n 

n 

t, , 

+ 

u ^ 

n t-^ 

+ •••  + u n 

j-1 

ij 

n2 

3-1 

j-1 

b.  Constraints  on  possible  combinations  of  processes  (6) 

(1-p  ) n t~]<i 

^ JEK* 

(1-P  )"^  n t < 1 .for  all  i 
JeK* 

c.  Constraints  on  possible  operating  ranges  (7) 

(l-Wy)  ty  i 1 

(1-PP^j)”^  t^j  £ 1 for  desired  i and  J 
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d.  Natural  constraints  (8) 

0 < t^j  ^ 1 for  all  1 and  j 


2.3  Model  Parameters.  As  indicated,  the  model  describes  a multistage 
activity,  the  stages  or  "reaches”  essentially  being  contiguous  sections  of  the 
stream.  The  activity  begins  at  a given  point  along  the  stream  at  a given 
time,  with  organic  wastes  being  first  treated  and  subsequently  deposited  into 
the  stream.  The  next  deposit  and  treatment  occurs  at  some  point  down  stream 
that  defines  the  end  of  the  given  stage  and  beginning  of  the  next. 


The  parameters  defining  the  model  are  mainly  those  Involved  in  the 
oxygen  profile  (sag)  equation.  For  a given  reach  at  a given  time  t the 
dissolved  oxygen  deficit  , given  in  the  previous  section,  is 

- Cj^L  + c^D  + Cj  (1) 

where  , 

c^  - K(r-K)’^e"‘^*'-e'*'‘') 

-rt 

C2  - e 

Cj  ■ (sludge  term) 


and  the  BOD  concentration  L of  the  stream  at  the  top  of  the  reach,  as 
indicated  in  the  last  section,  is 

V V 

T . ■ c t L. 

V’r  ^ «) 

and  BOD  concentration  at  the  bottom  of  Reach  1 is 

“ *=4  "-i  <3) 

-kiti 

where  c^  > e . The  parameters  Involved  in  those  equations  are  defined 

as  follows.  (For  notational  simplicity  the  subscript  i indicating  the 
reach  nuai>er  is  suppressed.) 
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K - Deoxygenation  constant.  In  day 
r - Reaeratlon  constant,  in  day  ^ 
t - Flow  time  along  the  reach,  in  days 
F - Fraction  of  river  bottom  covered  with  sludge 

Q - A coefficient  in  the  sludge  term,  determined  empirically 
in  day  ^ 

~ Oxygen  uptake  rate  per  unit  area  of  stream  bottom 
S 2 

surface,  in  gm/m  /day 

R - Hydraulic  radius  of  the  river  cross  section,  in 
meters 


V 


E 


E 


t 


e 


Volume  of  the  effluent  released  into  the  river, 
in  10^  gallons /day 

Concentration  of  effluent,  before  treatment, 
in  mg/f. 

Concentration  of  effluent,  after  treatment, 
in  mg/H 

Flow  rate  of  the  river  before  entering  the  reach, 
in  10^  gallons /day 


Aside  from  the  above  parameters  which  are  Involved  in  the  oxygen  sag 
equation  of  each  reach  the  following  parameters  will  also  enter  into  the 
formulation  of  the  model  at  various  stages. 

S - Maximum  allowable  oxygen  deficit,  in  mg/S. 

p - Minimum  required  or  feasible  BOD  removal  by  a specified 
sequence  of  treatment  components  in  a given  reach,  a 
fraction 

p - Maximum  possible  or  feasible  BOD  removal  by  a specified 
sequence  of  treatment  components  in  a given  reach,  a 
fraction 

PP  - Minimum  required  or  feasible  BOD  removal  by  a specified 
single  treatment  component  in  a given  reach,  a 
fraction 
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PF  - Maximum  possible  or  feasible  BOD  removal  by  a specified 
single  treatment  component  In  a given  reach,  a fraction 

1.  - Initial  BOD  level  of  the  stream  before  entering  Reach  1, 

0 In  mg/ 8, 

D.  - Initial  oxygen  deficit  of  the  stream  before  entering 
0 Reach  1,  In  mg/8 

3.  Formulation  of  the  Models  for 
the  Upper  Hudson  River 

The  pictorial  depiction  of  the  design  components  for  treatment 
facilities  and  corresponding  annual  operating  cost  data  considered  by  the 
author  of  the  model  along  the  Upper  Hudson  River  Is  given  In  Figure  4.  A 
description  of  the  treatment  components,  shown  In  the  figure  from  left  to 
right  Is  as  follows. 

Reach  1 - Primary  Clarifier  (PC) , Activated  Sludge  (AS) , 

Coagulatlon/Sedlmentatlon/Fllteratlon  (CSF) 

Reach  2 - PC,  Trickling  Filter  (TF),  Activated  Lagoon  (AL), 

CSF 


Reach 

3 

- PC, 

TF, 

AS, 

Carbon  Absorption  (CA) 

Reach 

4 

- PC, 

AL, 

CSF 

Reach 

5 

- PC, 

TF, 

AS, 

CSF 

Reach 

6 

- PC, 

TF, 

AL, 

CSF 

As  In  reference  [1],  under  three  different  treatment  policies,  three 
problems,  ♦ ^2  ’ ^3  respectively  formulated  which  Involve  the 

optimization  of  the  total  yearly  operating  cost  of  the  above  treatment  system. 

The  first  policy,  a fixed  dissolved  oxygen  requirement  policy,  requires 
that  a minimum  stream  dissolved  oxygen  concentration  level  of  6.2  mg/8  be 
maintained  along  the  entire  stream.  The  second  policy,  a uniform  BOD  treat- 
ment policy,  requires  that  each  treatment  plant,  (l.e.,  the  collection  of 
treatment  components  at  a given  discharge  point)  remove  at  least  95  percent 
of  the  BOD  concentration  of  Its  effluent.  The  third  policy  relaxes  the 
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f 

i 

i 


I 

I 

1 


\ 


I 

i 

* 

I 


I 

F 


requirement  of  the  first  two  policies,  but  uses  the  levels  of  dissolved 
oxygen  that  resulted  for  a given  reach  under  the  second  policy  as  the  minimum 
required  dissolved  oxygen  standard  for  that  reach.  This  allows  changes  in 
both  the  oxygen  deficit  and  treatment  levels  from  reach  to  reach,  while  guaran- 
teeing (with  the  data  as  given)  that  the  oxygen  standard  required  by  the  first 
policy  is  more  than  satisfied.  The  input  data  for  the  above  problems  were 
extracted  from  Reference  [1]  and  supplemented  in  part  by  J.  G.  Ecker,  the 
author  of  the  model.  The  formulations  of  Problems  P^,  P^*  and  P^  follow. 


Problem  P^:  Corresponding  to  a "fixed  dissolved  oxygen 

requirement  policy"  of  6.2  mg/£  along  the 
river. 

(Dimensions:  22  variables,  A2  constraints) 
Minimize  F - 19.4  + 86  t"'^®  + 152  t"*^^  + 

19.4  + 16.8  t"2*®®  + 27.4  t"'®^  + 179  t"^^^  + 

19.4  + 16.8  t32‘^®  + 91.5  t‘*^°  + 120  t"^®^  + 

19.4  + 45.9  tT;^^  + 179  tT;^^  + 

41  42  43 

19.4  + 16.8  + 91.5  t‘*^  152  t‘*^^  + 

19.4  + 16.8  + 27.4  t'l®^  + 179  t'*^^ 

61  62  63  64 


subject  to: 

a.  Dissolved  oxygen  deficit  constraint,  i.e., 

dissolved  oxygen  deficit  for  all  reaches  £ 2.18  mg/£  . 
(Type  (a).  Section  2.2) 


I 

I 
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1. 

'‘ll‘^ll*^12‘^ 

13^^ 

i 

1 

! 

r 

) 

2. 

“21 

+U22t 

21*^22 

*^23‘^2A-^ 

3. 

“31  " 

■"“32 

II 

^33' 

31*^32‘^33‘^3A^^ 

1 

A. 

“A1  " 

■^A2 

It 

■^A3 

■^AA*^Ai‘^A2 

^3-^  ; 

5. 

“51  " 

+U52 

II 

•^53 

" ■^“5A  " 

■^55*^51*^52‘^53*^5A- ^ j 

6. 

“61 

■^62 

II 

■^63 

" -^bA  " 

■^65  " ■*^66‘^61*^62*^63^6A-^  | 

b.  Constraints  on  combination  of  processes 
(Type(b),  Section  2.2) 


1 

] 

7. 

.1 

h2  1 ^ 

at 

on 

1 

i 

8. 

.15 

‘ai 

^22  *^23  - ^ 

at 

on 

1 

j 

9. 

.15 

*^32  *^33  - ^ 

at 

on 

\ 

] 

10. 

.15 

^aU^ 

at 

on 

i 

i 

1 

11. 

.15 

si 

S2  ^53  ^ ^ 

at 

on 

t 

i 

i 

12. 

.15 

si 

*^62  *^63  - ^ 

at 

on 

:! 

• \ 

13. 

1.A275  t 

61  ^62  1 1 

at 

In 

on  reach  6 


c.  Operating  range  constraint  for  components 

of  treatment  facilities  (Type  (c).  Section  2.2) 


lA.  1.25  t,,  < 1 
A1  — 


at  least  20  percent  removal  by  first  component 
In  reach  A 
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d.  Redundant  constraints  (added  to  prevent 
numerical  overflow) 


-80 

-1 

-1 

-1 

15. 

e 

hi 

h2 

*^13 

< 

1 

-80 

-1 

^-1 

^-1 

^-1 

16. 

e 

*^21 

*^22 

*^23 

24 

< 

1 

-80 

-1 

-1 

-1 

-1 

17. 

e 

hi 

*^32 

*^33 

*^34 

< 

1 

-80 

-1 

-1 

-1 

18 . 

e 

*^41 

h2 

*^43 

< 

1 

-80 

-1 

^-1 

^-1 

^-1 

19. 

e 

hi 

*^52 

h3 

•^54 

< 

1 

-80 

-1 

^-1 

-1 

-1 

20. 

e 

*^21 

*^22 

ft 

N> 

1 

h4 

< 

1 

e.  Natural  constraints  (Type  (d) , Section  2.2) 

21-42.  0 < t^j  < 1 1-1,6,  j - 1 , 

The  coefficients  u^^j  In  the  oxygen  deficit  constraints  (1-6)  are 

complicated  and  lengthy  functions  of  nearly  all  the  problem  parameters.  To 
conduct  a detailed  sensitivity  analysis  (the  subject  of  a follow-up  report) 
with  respect  to  the  problem  parzuneters  It  was  necessary  to  delineate  the 
functions  defining  these  coefficients  In  terms  of  the  problem  parameters. 
Explicit  derivation  of  these  functions  by  hand  proved  extremely  tedious  and 
time  consuming  and  evaluations  were  prone  to  computational  errors.  The  task 
was  virtually  prohibitive  for  more  than  three  reaches,  since  the  equations 
for  a given  reach  Involve  quantities  that  depend  on  results  obtained  In  all 
previous  reaches.  To  resolve  this  difficulty,  the  computer  program  listed 
In  Appendix  2 was  written  to  define  and  calculate  the  u^j  coefficients  In 

terms  of  the  problem  parameters.  The  computed  u^j  for  Problem  P^,  along 

with  a listing  of  the  parameter  values  are  given  in  Appendix  2.  Although  the 
program  was  written  to  handle  the  problem  Involving  six  reaches  and  the 


Upper  Hudson  River  data  base,  It  can  readily  be  modified  to  handle  other 
data  and  other  likely  variations  that  might  appear  in  other  models. 


Problem  P^:  Corresponding  to  a "uniform  BOD  treatment 

policy"  of  95  percent  removal  for  all  reaches. 
(Dimensions:  22  variables,  48  constraints) 


Minimize  F (See  Problem  P^) 


Subject  to: 


a.  Constraints  1-20  of  Problem  Pj^ 

b.  At  least  95  percent  treatment  In  all 
reaches,  l.e.. 


21.  20  < 1 


22.  20  < 1 


23.  20  tjj  tjj  tjj  < 1 


“•  “ '41  '42  '43  i ' 


25.  20  tjj  tjj  tjj  < 1 


'61  '62  '63  '64  i ' 


c.  Monnegatlvlty  constraints  of  Problem  P^ 


Based  on  the  assumptions  and  results  of  this  problem,  the  dissolved 
oxygen  level  (hence  the  oxygen  deficit  S^)  in  each  reach  was  calculated. 

The  oxygen  deficit  levels  so  obtained  were  used  as  the  new  oxygen  deficit 
standards  for  Problem  P^  and  new  u^  coefficients  were  calculated.  Thus, 
Problem  P^  was  formulated  as  follows. 
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Problem  P^:  Corresponding  to  the  policy  of  using  the  levels 

of  dissolved  oxygen  that  resulted  for  a given 
reach  for  Problem  ^2  minimum  required 

dissolved  oxygen  standard  for  that  reach. 
(Dimensions:  22  variables,  42  constraints) 


Minimize  F (See  Problem  Pj^) 

Subject  to: 

a.  Constraints  1-6  of  Problem  Pj^  with  coefficients 
u^j  calculated  using  new  S^^ 


b.  Constraints  7-20  of  Problem  P^^ 


c.  Nonnegativity  constraints. 


The  coefficients  u^^  along  with  the  required  data  for  this  calculation 
are  listed  in  Appendix  4. 

By  the  usual  transformation,  l.e. , letting 

6 

X.  » -f,n  t , (i-1,6,  j"l,m.  and  k*l,  E m ) , 

as  labelled  in  Figure  4,  the  geometric  programming  problems  » ^2  ’ 

P^  were  converted  to  convex  progranmlng  problems.  The  resulting  equivalent 
problems,  labelled  ^^,02*  and  C^  , respectively,  follow. 


Problem  C^:  Convex  Equivalents  of  Problem  Pj^. 

(Dimensions:  22  variables,  42  constraints) 


1.47x, 


Minimize  F - 19. 4e 


+ 86e 


. 38x-  .27x. 

^ + 152e 


1.47x, 


+ 19. 4e 


+ 16. 8e 


1.66x,  •37x_ 

^ + 27. 4e  **  + 179e  ' 
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1.47xg  1.66Xq  .33x. 

+ 19. 4e  + 1.68e  ^ + 91. 5e  ^ + 120e 

1.47x,„ 

+ 19.4e  + 45.9e  + 179e 


1.47x^ 

+ 19. 4e  16. 8e 

1.47x^ 

+ 19. 4e  16. 8e 


1.66x, - •30x,t  .27x, „ 

91. 5e  + 152e 


1.66x-^  •^^*59 

27. 4e  + 179e 


subject  to: 


a.  Dissolved  oxygen  deficit  constraint,  l.e.,  the 

dissolved  oxygen  deficit  for  all  reaches  £ 2.18  rag/ 2 . 
(Type  (a).  Section  2.2) 


1.  + 1 

2.  -«2j 

3.  -«3,  " 

-32 

4.  ^,3  - 

-u  " 

42 

-43 

:<»12*hs^l4> 
-44*  ♦ 

5.  -53  - 

-n  " 

—52 

-53 

—54  " — 55‘ 

4.  -^3  - 

— ” 

— *2 

-43 

—44  " —45 

-(*,,+*  ,+«  4*  ) 


i 


b.  Constraints  on  combination  of  processes 
(Type  (b).  Section  2.2) 

~(xj^'^2^  + AnlO  ^ 0 (At  most  90  percent  reoioval  by  total  secondary 

treatments  in  Reach  1.) 
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-(x,+x-+x,)  + £.n  TT  > 0 
456  .15  — 


-(x8+X9+x^o)+  «-«  .IJIO 


-(X12+X13)  + .l5  > 0 


-(X15+X18+X12)  + Jin  > 0 


-(Xi9+X2o+*2i)  + .l5  ^ ° 


(At  most  85  percent  removal  by 
total  secondary  treatments  in 
Reaches  2 - 6 . ) 


13.  * Q + *00  + -7  > 0 (At  least  30  percent  removal  by  first  two 

components  In  Reach  6.) 

c.  Operating  range  constraint  for  components 
of  treatment  facilities 

(Type  (c) , Section  2.2) 

14.  X,,  + Jin  .8  0 (At  least  20  percent  removal  by  first 

component  in  Reach  4.) 

d.  Redundant  constraints  added  to  prevent 
numerical  overflow. 


15. 

-(Xj^+X2+X2) 

+ 80  > 0 

16. 

-(x^+x^+x^+x^) 

+ 80  > 0 

17. 

-(Xg+x^-hCio+Xii) 

+ 80  > 0 

18. 

-(Xi2+Xi3+Xi^) 

+ 80  > 0 

19. 

+ 80  ^ 0 

20. 

-(Xi9-fx2o+*2l‘^22>  + 80  > 0 

t 

e. 

Natural  constraints 
(Type  (d).  Section  2.2) 

1 

> 0 i 

- 1,22 
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Problem  C^: 

Minimize  F 
subject  to: 

a. 

b. 


Convex  Equivalent  of  Problem  P^ 

(Dimensions:  22  variables,  48  constraints) 

(See  Problem  C^) 

Constraints  1-20  of  Problem  C^ 

At  least  95  percent  treatment  in  all 
reaches,  l.e.  , 


Xi  + X2  + X3  + S-n  .05  ^ 0 

X,  + X-  + X,  + x_  + in  .05  > 0 

4 5 o 7 — 

^8  *9  ^ *10  ^ *11  ^ ® 

Xj^2  *j^3  *34  ® 

. + x.,  + X- T + x- - + in  .05  > 0 

> 16  17  18 

, + ^20  *21  ^ *22  ^ ° 


c.  Nonnegativity  constraints 


Problem  C^  : 


Maximize  F 
subject  to: 


b. 


Convex  Equivalent  of  Problem  P^ 

(Dimensions:  22  variables,  42  constraints) 

(See  Problem  C^) 

Constraints  1-6  of  Problem  C^  with 
coefficients  u^^j  of  Problem  P^ 

Constraints  7-20  of  Problem  C^^ 


c.  Nonnegativity  constraints 
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4.  Optimal  Solutions  of  the  Upper 
Hudson  River  Models 

The  convex  prograomlng  problems  C^.  and  developed  In  the 

previous  section  were  solved  with  SUMT  for  the  optimum  treatment  levels 
and  Che  annual  waste  treatment  cost.  The  approximate  resulting  optimal 
solutions  obtained  for  the  respectively  equivalent  problems,  P^,  and 
are  as  follows. 

4.1  Problem  Pj^ 

• Optimal  treatment  levels  (using  a fixed  dissolved 
oxygen  standard  of  6.2  mg/1  along  the  entire  stream) 


t* 

1-1 

j-1 

.7236 

j-2 

.4095 

j-3 

.6617 

j-4 

1-2 

.7729 

.7854 

.2471 

.4343 

1-3 

.8640 

.8668 

.4303 

1. 

1-4 

.8000 

.3229 

1. 

— 

1-5 

1. 

1. 

1. 

1. 

1-6 

.7729 

.7854 

.2471 

.9551 

* Annual  waste  treatment  cost  > $1,837,851 

• Binding  constraints:  (6,8,12  and  14) 

The  fact  that  a given  component  t^^  is  equal  to  unity  at  the 

optimal  solution  Indicates  that  the  corresponding  treatment  facility  need 
not  be  included  (if  this  Is  operationally  feasible)  In  the  plant  design, 
since  that  component  contributes  nothing  to  the  treatment  process  when  It 
Is  operated  optimally.  The  calculated  annual  waste  treatment  cost  should 
then  be  reduced  accordingly.  In  the  above  table,  these  components  were 
not  removed  in  calculating  the  total  costs.  Suppose,  however,  that  we 
conclude  that  the  last  treatment  facilities  In  Reaches  3 and  4 and  all 
facilities  in  Beach  5 are  not  needed.  Making  this  adjustment,  the  optimal 
waste  treatment  cost  as  shown  In  Table  1 Is  $1,253,140  per  year. 
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Table  1 


OPTIMAL  PLANT  BOD  REMOVAL  LEVELS 
AND  COSTS  FOR  PROBLQ!  P^^ 


(Using  a fixed  dissolved  oxygen  standard 
of  6.2  mg/S.  along  the  entire  stream) 


Reach 

Mo. 

Total  BOD 
Removal 

(X) 

Cost 

$/Year 

Redundant 

Components 

1 

80.39 

321,870 

None 

2 

93.48 

363,230 

None 

3 

67.77 

163.190 

CA 

4 

74.17 

103,260 

CSF 

5 

0.00 

— 

All  Components 

6 

85.67 

301,590 

None 

Total 

$1,253,140 

If  we  exclude  the  last  component  In  Reach  6,  as  Is  done  In  [1],  the 
cost  would  further  decrease  to  $1,253,140  - $182,060  ■ $1,071,070  per  year. 
(The  solution  obtained  In  [1]  Whs  $1,106,000  per  year.)  However,  based  on 
our  results  we  do  not  see  the  Justification  for  the  given  exclusion. 

As  mentioned  previously,  we  have  traced  back  the  discrepancy  of  our 
solution  with  that  obtained  in  [1]  to  the  differences  In  the  values  of  the 
coefficients  u^^  used  to  calculate  the  oxygen  deficit  constraints.  Our 

calculated  u^j  , along  with  the  list  of  required  Input  data,  are  Included 

In  the  Appendix. 


optimal  treatment  levels  (using  95%  uniform 
treatment  for  all  reaches) 


t* 

111 

i-1 

J-1 

.6338 

j-2 

.2452 

j-3 

.3216 

J-4 

i-2 

.7729 

.7854 

.2471 

.3333 

i-3 

. 7720 

.7845 

.2477 

.3333 

i-4 

.7792 

.2057 

.3333 

— 

1-5 

.7720 

.7845 

.2477 

.3333 

1-6 

.7730 

.7854 

.2471 

.3333 

Optimal  annual  waste  treatment  cost  » $2,330,279 

Binding  constraints:  (8,9,10,11,12,21,22,23,24, 
25,26) 


Unlike  the  results  obtained  for  Problem  P^,  none  of  the  components 
of  the  solution  matrix  are  unity.  Thus,  we  cannot  exclude  any  treatment 
component  with  this  policy.  As  shown  in  Table  2,  if  we  exclude  the  treat- 
ment facilities  of  Reach  5,  as  is  done  in  Reference  [1],  then  the  optimal 
annual  waste  treatment  cost  would  reduce  to  $2,330,279  - $397,080  ■ $1,933,199, 
which  is  in  close  agreement  with  the  result,  $1,926,000  per  year,  given  in 
Reference  [1].  (In  fact,  if  we  use  our  Table  2 calculation,  we  get 
$2,322,580  - $397,080  * $1,925,900,  which  agrees  almost  exactly  with  [1].) 
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Table  2 

OPTIMAL  PLANT  BOD  REMOVAL  LEVELS 
AND  COSTS  FOR  PROBLEM  P^ 

(Using  a 95%  uniform  treatment 
for  each  reach) 


Reach 

No. 

Total  BOD 
Removal 
(%) 

Cost 

$/Year 

Redundant 

Components 

1 

95 

391,120 

None 

2 

95 

388,300 

II 

3 

95 

365,030 

II 

4 

95 

393,150 

II 

5 

95 

397,080 

11 

6 

95 

388,300 

It 

Total  $2 

,322,980* 

*Thls  value  Is  about  $7,300  less  than  that  found  by  direct  solution.  This 
discrepancy  Is  attributed  to  the  translation  errors  Involved  In  hand  calculation 
of  the  costs  on  the  above  table. 

All  of  the  uniform  treatment  constraints,  l.e.,  constraints  21-26,  are 
binding.  On  the  other  hand,  the  dissolved  oxygen  constraints  1-6  are  more 
than  satisfied.  As  In  [1],  this  suggests  using  the  dissolved  oxygen  levels 
resulting  under  the  uniform  treatment  policy  as  the  minimum  standard  for  each 
reach,  thus,  allowing  the  levels  of  treatment  to  vary  from  reach  to  reach. 

Under  such  conditions  one  would  hope  that  the  optimal  cost  would  be  lower  than 
that  obtained  under  the  policy  of  uniform  treatment  while  yielding  an  optimal 
treatment  schedule  that  would  be  nearly  uniform.  (Uniform  policies  are 
apparently  easier  to  administer.)  The  following  optimal  solution  of  problem 
P^  confirms  these  expectations. 
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As  it  Is  seen,  the  last  component  in  reach  5 is  redundant  since 

t-  . « 1 . Thus  the  optimal  waste  treatment  cost 
5,4 

$2,147,336  - $152,000  ■ $1,995,336  per  year.  The  solution  obtained  in 
Reference  [1]  was  $1,460,000  per  year.  Here  again,  the  discrepancy  between 
our  solution  and  that  of  Reference  [1]  is  partly  due  to  the  differences 
in  the  values  used  as  the  coefficients  u^j  of  the  dissolved  oxygen  con- 
straints, and  partly  due  to  the  fact  that,  for  reasons  explained  in  Section 
4,1,  we  have  not  excluded  the  treatment  facilities  in  Reach  5 in  our 
calculations. 
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Table  3 

OPTIMAL  PLANT  BOD  REMOVAL  LEVELS 
AND  COSTS  FOR  PROBLEM  P^ 


Reach 

No. 

Total  BOD 
Removal 

m 

Cost 

$/Year 

Redundant 

Components 

1 

95.02 

391,300 

None 

2 

97.33 

458,380 

II 

3 

88.40 

323,210 

It 

4 

86.02 

308,090 

II 

5 

21.04 

134,390 

CSF 

6 

94.12 

372,650 

None 

Total  $1,988,020* 

*This  value  is  about  $7,j16  less  than  that  found  by  direct 
solution.  This  discrepancy  is  attributed  to  the  truncation 
errors  involved  in  hand  calculation  of  the  costs  in  the  above 
table. 

As  it  is  seen  in  Table  3,  the  treatment  levels  for  all  reaches 
except  Reach  5,  which  contributes  a moderate  level  of  pollutant  to  the 
stream,  are  nearly  uniform.  Also  the  optimal  treatment  cost  is  $1,995,336 
per  year  compared  to  that  for  uniform  treatment  of  $2,330,279  per  year. 
Therefore,  under  the  stream  standards  (dictated  by  the  solution  of  P2) 
we  may  achieve  a relatively  uniform  treatment  policy  and  at  the  same  time 
save  $2,330,279  - $1,995,336  - $334,943  per  year. 

A comparison  of  the  above  results  shows  that  the  minimum  annual  waste 
treatment  cost  is  obtainable  under  the  fixed  dissolved  oxygen  policy  of  6.2 
mg/Jl  . The  cost  incurred  under  the  95  percent  uniform  treatment  policy  is 
by  far  the  highest,  although  this  policy  is  the  easiest  to  monitor  and  ad- 
minister. The  same  stream  standards  achieved  under  uniform  treatment  level 
policy  can  be  met  under  an  almost  uniform  treatment  level  policy,  with  appre- 
ciable savings  in  treatment  cost.  The  annual  saving  under  such  a policy 
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for  the  model  under  study  is  of  the  order  of  $335,000,  amounting  to  a 14.4 
P'‘'‘cent  savings  In  annual  treatment  cost. 

5.  Conclusions 

Based  on  the  previous  results  and  discussions,  the  following  conclusions 
are  Inmediate. 

* Among  the  three  different  policies  studied  for  optimal 
water  pollution  control  and  treatment  facility  allocation, 
the  policy  requiring  a fixed  dissolved  oxygen  of  6.2  mg/Jl 
along  the  stream  yields  the  least  annual  waste  treatment 
cost,  $1,253,140. 

* The  requirement  of  a 95  percent  uniform  waste  treatment 
level  In  all  reaches,  a policy  easier  to  monitor  and 

administer,  requires  an  additional  annual  coat  of  $1,077,140  ! 

I 

when  this  requirement  Is  added  to  the  fixed  oxygen  level  i 

standard. 

* The  same  dissolved  oxygen  standards  achieved  under  95  percent 
uniform  treatment  level  policy  can  be  met  under  a near-uniform 
treatment  level  policy,  with  savings  of  about  $335,000  per 
year,  amounting  to  a 14.4  percent  saving  In  the  annual  treat- 
ment cost  obtained  by  a uniform  treatment  policy. 

* The  above  findings  are  consistent  with,  but  not  identical  to, 
the  findings  reported  in  [1].  The  reason  for  discrepancies 
in  the  optimal  solution  is  attributed  to  the  discrepancy 

between  u , the  coefficients  appearing  in  the  dissolved  | 

oxygen  constraints,  used  by  these  two  studies.  1 

* The  computer  routine  SUMT  solved  these  geometric  prograimning 
convex  equivalent  models  without  difficulty.  The  efficiency, 
however,  compared  to  that  of  a special  purpose  GP  code  for 
solving  the  original  problem,  has  yet  to  be  determined.  Apart 
from  this,  the  given  model  can  now  be  appreciably  modified 
to  Include  much  more  general  functions  (than  those  required 
by  the  G.  P.  model),  since  SUMT  assumes  no  such  special 


structure. 
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APPRNDIX  1 


Dissolved  Oxygen  Deficit  Constraint  for  Reach  1 In  Problem  P^^  . 

The  oxygen  deficit  D,  at  the  bottom  of  Reach  1 Is  calculated  from 

'’l 

Equation  (1)  (Section  2.3)  as  follows 

°b^  " ^11  4 ^21  °bQ  ^31 


where 


Cii  - Ki(vKp-l|.  'I'l) 


Si-' 


“’^1*^1 


Si  ■ SS  sji-' 


and 


* V — K \(using  Equation  2,  Section  2.3) 

\ *1  % '’o/ 


s ■ S J-  Sj  • 

j“l 

Imposing  the  dissolved  oxygen  deficit  constraint  at  the  bottom  of 
Reach  1,  i.e.,  D,  < S,  , yields 

bi  - 1 


^11 


\*\  ”"0  ^ '''  - ''  ■ 


After  simplification,  we  obtain 

V, 


E, 


R,  W S -C,,b„-C,  -C, 


11  +V„  il  I 1 21“0  31  11  V„  +V 


"l  ‘'o 


^1  V 


ij  - 
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using  the  data  provided  In  Appendix  3.  The  coefficient  of  H t.. 

j-1  ^ 

In  the  above  constraint  (which  Is  labelled  In  the  model)  was 

calculated  to  be  0.7756.  Thus,  the  oxygen  deficit  constraint  for  Reach  1 
becomes  simply 

.7756  < 1 


Dissolved  Oxygen  Deficit  Constraint  for  Reach  2 In  Problem  . 

The  oxygen  deficit  D,  at  the  bottom  of  Reach  2 Is  calculated 
from  Equation  1 (Section  2.3)  as  follows: 


where 


^2  ’ "l2  "*2  ^ "22  °b,  + "32 


"22-® 


-r_t- 
2 2 


"32  - Q2"2‘'s. 


'2  V. 


"I"?; — Vv  ^ \(“8lng  Equation 


2,  Section  2.3) 


®2  - h *^23 


'^R  - '^R  ''e 

1 *^0  ®1 


-K^t 

- e (using  Equation  3,  Section  2.3) 

L,  ■ i,: — (v„  E,  n t..  + V_  L.  \ 

^ ^0  ®oy 
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INPUT  DATA  AND  CALCULATED 
FOR  PROBLEMS  Pj^  AND  P2 
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Appendix  4 

INPUT  DATA  AND  CALCULATED 
FOR  PROBLEM  P^ 


THE  GEORGE  WASHINGTON  UNIVERSITY 
Program  in  Logistics 
Distribution  List  for  Technical  Papers 


Th«  C«ors«  WMhinftoA  Unlvcratcy 
0((ic«  of  SpoiMorod  ftooMrch 
Library 

Vlca  Praal4ant  H.  f.  Bright 
Da«a  Harold  Llabowlti 
Da«n  Hanry  Soloaoo 

OMR 

Chiaf  of  Naval  laaaarch 
(Codaa  200,  41A) 

Raaidaot  Mpraaantattva 

OPHAV 

OP-40 

DCNO,  Logiatlca 
Navy  Dapt  Library 
NAVDATA  AutoaMtlon  Cad 
0P-9b4 

Naval  Aviation  Intagratad  Log  Support 

NAADAC  Tach  Library 

Naval  Clactronlca  Lab  Library 

Naval  Pacllltlaa  Eng  Cad  Tach  Library 

Naval  Ordnanca  Station 
Loulavllla,  Ky. 

Indian  Haad,  Hd. 

Naval  Ordnance  Sya  Cad  Library 

Naval  Raaaarch  Branch  Of flea 
Bouton 
Chicago 
Naw  York 

Paaadana 

San  rranclaco 

Naval  Ship  Eng  Cantar 
Philadelphia.  Pa. 

Uaahlngton.  DC 

Naval  Ship  fUa  4 Dav  Cantar 

Naval  Saa  Syataaa  Conaand 

PKS  30611 
Tach  Library 
Code  073 

Naval  Supply  Syaeaa«  Coaaand 
Library 

Oparaclona  and  Inventory  Analyjle 

Naval  War  Collage  Library 
Navport 

BUPEIIS  Tach  Library 
PMSO 

Integrated  Saa  Lift  Study 

USN  kmo  Depot  Earle 

USN  Poatgrad  School  Nontarey 
Library 

Dr  Jack  R.  Boratlng 
Prof  C.  R.  Jonaa 

US  Narine  Corpa 
Co— andant 

Deputy  Chiaf  of  Staff*  R4D 

NarlM  Corpa  School  Quant  Ico 
Landing  Perea  Dav  Ctr 
Logiatiea  Officer 
Co— ending  Officer 

USS  Prancla  Marion  (LPA-249) 

Ar— d Percaa  Induatrtal  Collage 

Ar— d forcaa  Staff  CoJlag* 

Ar«y  War  Collage  Library 
Carliala  Barraelta 

Aray  Cod  4 Can  Staff  Collaga 


Arny  Logiatlcv  Hgt  Canter 
Port  Lee 

Co— ending  Officer.  USALOSKA 
New  Cuobcrland  Aroy  Depot 

Arny  Inventory  Rea  Ofc 
Phllade Lphla 

Air  Force  Headquertera 
AFADS> 1 
LEXY 
SAF/ALC 

Griff iae  Air  Force  Base 

Rallabillty  Analyala  Cantar 

Gunter  Atr  Force  Base 
AFUMC/XR 

Maxwell  Air  Force  Base  Library 

Wrlght'Pattereon  Air  Furcc  Boat- 
Log  ConoMnd 
Reaearch  Sch  Log 
AFALD/XR 

Dcfenae  Docunentat Ion  Center 

National  Acadeny  of  Science* 

Haritlne  Trananortation  Kee  Board  Library 

National  Bureau  of  Standards 
Dr  B.  H.  Colvin 
Dr  Joan  Koaenblatt 

National  Science  Foundation 

Nailonal  Security  Agency 

Weapon  Syrtent  Evaluation  Group 

brttUh  Navy  Staff 

National  Uefenae  Hdqtrs,  Ottawa 

LoKlstiCK.  OR  Analyala  hHtabllahncnt 

American  Power  Jet  Co 
George  Chcrnowltt 

General  Dynaalcs,  Pomona 

General  Reaearch  Corp 
Or  Hugh  Cole 
Library 

Loglactca  Management  Inetltute 
Dr  Murr.iy  A.  Geieler 

MATHTEC 

Dr  Eliot  Feldman 

Hand  Corporal  Ion 
Library 

CarneglcoHellon  University 
Dean  H.  A.  Simon 
Prof  c.  Thoaipaon 

Caae  Western  Reserve  University 
Prof  B.  V.  Dean 
Prof  N.  Mesarovlc 
Prof  S.  Racks 

Cornell  University 

Prof  K.  E.  Bechhofer 
Prof  R.  Conway 
Prof  Andrew  Schulte.  Jr. 

Cowles  Foundation  for  Research  in  Ecen— ice 
Prof  Herbert  Scarf 
Prof  Martin  Shublk 

Florida  State  University 
Prof  It.  A.  Bradley 

Harvard  u. .veralty 
Prof  K,  J.  Arrow 
Prof  W.  C.  Cochran 
Prof  Arthur  Schtelfer.  Jr. 

Princeton  University 
Prof  A.  W.  Tucker 
Prof  J.  w.  Tukey 
Prof  Geoffrey  S.  Watmon 


Purdue  University 
Prof  S.  S.  Cupte 
Prof  H.  Rubin 
Prof  Andrew  Whlnecon 

Stenford  University 
Prof  T.  U.  Anderson 
Prof  C.  B.  Dentslg 
Prof  P.  S.  Hllller 
Prof  0,  L.  Iflehert 
Prof  SeMiel  Kerlln 
Prof  G.  J.  Llebernan 
Prof  Herbert  Solonon 
Prof  A.  P.  Velnott,  Jr. 

University  of  Cellfornls.  Berkeley 
Prof  R.  B.  Berlov 
Prof  0.  Gele 
Prof  Jeck  Kiefer 
Prof  Rosedith  Slcgreeves 

University  of  Celifomle,  Los  Angeles 
Prof  J.  R.  Jeckson 
Prof  R.  R.  O'Neill 

University  of  North  Csrollne 
Prof  U.  L.  Snlth 
Prof  H.  R.  Leadbetter 

University  of  Pennsylvania 
Prof  Russell  Ackoff 
Prof  Thonas  L.  Saaty 

University  of  Texas 
Prof  A.  Charnea 

Yale  University 

Prof  P*  J.  Anaconbe 
Prof  I.  R.  Savage 

Prof  Z.  W.  Blrnbaun 

University  of  Washington 

Prof  B.  H.  Blaslnger 

Tlie  Pennsylvania  State  UniversUy 

Prof  Seth  Bonder 

University  of  Michigan 

Prof  G.  E.  P.  Box 

University  of  Wisconsin 

Or  Jeroflie  Bracken 

foscftute  for  ■'iefeitae  Anmlyn^t 

Prof  H.  Chernoff 

Hats.  Institute  of  Technology 

Prof  Arthur  Cohen 

Rutgers  > The  State  University 

Mr  Wallace  M.  Cohen 

US  General  Accounting  Office 

Prof  C.  Dentan 

Coluebla  University 

Prof  Hasao  Pukushlaia 
Kyoto  University 

Prof  Saul  I.  Casa 

University  of  Maryland 

Dr  Donald  P.  Caver 
Carnel*  California 

Prof  Aairit  L.  Coal 
Syracuse  University 

Prof  J.  F.  Hannan 

Michigan  Stats  University 

Prof  H.  0.  Hartley 

Texas  ABM  Foundation 


Mr  Gerald  F.  Hein 

NASA,  Lewis  Research  Center 

Prof  W.  M.  Hirsch 
Courent  institute 

Dr  Alan  Hoffaan 

IBM,  Yorktown  Heights 

Prof  John  R.  Jsbell 

State  University  of  New  York,  Aad^erst 

Dr  J.  L.  Jain 

University  of  Delhi 

Prof  J.  H.  K.  Kao 

Polytech  Inatitute  of  New  York 

Prof  W.  Kruakal 

University  of  Chicago 

Hr  S.  KuaMT 

University  of  Hedrss  . 

Prof  C.  E.  Leaks 

Rensseleer  Polytsch  Inetituce 

Prof  Loynes 

University  of  Shefflsld,  England 

Prof  Steven  Nehaiiee 

University  of  Pittsburgh 

Prof  0.  B.  Owen 

Southern  Hethodiat  Uolverelty 

Prof  E.  Person 

Texas  ABM  University 

Prof  H.  0.  Postsn 

University  of  Connecticut 

Prot  R.  Renege,  Jr. 

University  of  Delewers 

Prof  Hans  Rledwyl 
University  of  Bern 

Dr  Pred  Rigby 

Texas  Tach  Collage 

Hr  David  Rosenblatt 
Weehlngton,  D.  C. 

Prof  H.  Roeenblett 

Unlverolty  of  Celifomle,  Sea  Diego 

Prof  Alan  J«  Rowe 

Unlveraity  of  Southern  Celifomle 

Prof  A.  H.  Rubenstein 

Northwestern  University 

Or  H.  E.  Selveoon 
West  Loa  Angelea 

Prof  Edward  A.  Sliver 

Univerelty  of  Waterloo,  Ceaede 

Prof  N.  J.  Sobel 

Georgia  last  of  Techaelog;y 

Prof  I.  H.  Thrall 
Rice  Univerelty 

Dr  S.  Vejda 

Uelvereity  of  Sussex,  tatlead 

Prof  T.  H.  Wbltin 
Ueeleyen  University 

Prof  Jacob  Wolfowlta 

Univerelty  of  South  Ploride 

Prof  Hex  A.  Woodbury 
Dufco  Univeroity 


Dece^r  tSTI 
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FNC/NflklNG  HOPIS  M)P  IHl  loMUUPONA'S  AS  vv'plMIN  IN  IMl.  Kl.ti»Wl>S  1)1 
FOlLOWIsr,  /-.Ol  I HNMl  M AL  AND  1>  H 1)  t h S S 1 U N A 1.  INClNll.WlNl.  DM  C,  A N V / A I IDNS  AMU 
THOSf  OI  I H I S f;  L OR  C r WASH  I FiCT  l)N  II  N 1 V I W.*.!  I Y 

IHliUU  m t OMMISSIONhHS  I)  I S I M K I 01  CDLISIOIA  • i 

U M T I 1)  SI  Ills  AlOMIC  r N F R rOMMlSSION 

DlPARIsifNI  Ol  IMF  ARMS  lINflKl)  SVATf.S  01  AMlltltA  . ^ 

DtPAHI'IINl  OF  IHt  NAVY  IlNIIFD  S 1 A 1 F.  S O I AMLRlCA 
DFPAkl'IF.SI  OF  IHE  AIR  FORCE  II N 1 f f.  U S:  1 A 1 1 S 01  AMI  RICA 

NAiioNiL  advisory  cosimiftiv:  Fob  aironai  ncs 

NAIIONAi.  BintAV  OF  SJaNDaRDS  11  S DFPAHIMINI  OI  aos'mIRCI 

sorii  IS  of  civii  i ncinf.ibs 

A All  RICAN  INSIIII  It.  UF  I 1,  F'(  I R I C AL  r NG  I Ml  HS 
FHF  A Sit  RICAN  StlCll.rV  OF  MECHANICAL  1 N C,  I M E R S 
IHF  SOCIF'rV  OF  AMFFiICAN  MILIIARS  1 N C,  1 N I I R S 

AMI  RICAN  INSIIII  IE  OF  MINING  fV  M b.  1 A L L 11  R C.  I C A L ENCIM.IRS 
' DISIRU  I Ol  COILMIIIA  SOCIfirs  Of  PROFESSIONAL  1.  N G I M E R S,  I \C  ^ ' 

» t li  I I S S n I L I I OF  R A DIO  1 NG  I N F,  I R S I NC  '«*?■' 

riFF  CHFMICaL  INGINEERS  C L U li  Ol  WASHINCVON 

washingfon  son  e is  oi  enginifhs  i 

r fVIlKNFR  KINGSHLRV  a S I F.N  FloUSK  - A WC  H 11  1 C 1 S ■ Kt-. 

CHAfiEF  S H lOSfPKlNS  COM  P AN  ■>  - mi  I I DF  liS  , , ' J 

.SFbCFETS  OF  WOSIFN  ENGINEERS  ^ ' , *• 

' Nil  ION  AE  ACADFMV  of  SCIENCES  NAIIONAE  R I St  ARCH  ( O i S f 1 L«  v^' 

rUF  PURPOSE  OF  THIS  VAULT  IS  INSPIRTTJ  ^ Df.lMCWlV)  TO 

CHARLES  HOOK  TO'MPMNSv  iioTjou  or  i NG  I s b r wfnc. 
ojflfCAr.sf  OF. ms  f n r,  i n ee:r  i nc  con! ’r  i lo  O.o's'jS  r i!i  s . f's'i T R sr  * y i o f'ms 


>~'Jrhsi  \tl  .MJXfjfo  HIS  NAIION  and  lo  ^Th  r^<*  1 s 

^ BY  ThE  GEORGE  W'ASHINCYoN  IINIVE.KSITV 

•■'ROHEftT  a'-^ElikMlNC  ‘ fx  Cl  OV:>  H-  sr  \HV\N 


VI  BO  A fc  H o » i ij  I y n f ^ 


To  cope  with  the  expanding  technology,  our  society  must 
be  assured  of  a continuing  supply  of  rigorously  trained 
and  educated  engineers.  The  School  of  Engineering  and 
Applied  Science  is  completely  committed  to  this  ob- 
jective. 


